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We report that the tt band of graphene sensitively changes as a function of an external potential 
induced by Na especially when the potential becomes periodic at low temperature. We have mea- 
sured the band structures from the graphene layers formed on the 6H-SiC(0001) substrate using 
angle-resolved photoemission spectroscopy with synchrotron photons. With increasing Na dose, the 
TT band appears to be quickly diffused into background at 85 K whereas it becomes significantly en- 
hanced its spectral intensity at room temperature (RT). A new parabolic band centered at k ^1.15 
also forms near Fermi energy with Na at 85 K while no such a band observed at RT. Such 
changes in the band structure are found to be reversible with temperature. Analysis based on our 
first principles calculations suggests that the changes of the tt band of graphene be mainly driven by 
the Na-induced potential especially at low temperature where the potential becomes periodic due 
to the crystallized Na over layer. The new parabolic band turns to be the tt band of the underlying 
buffer layer partially filled by the charge transfer from Na adatoms. The five orders of magnitude 
increased hopping rate of Na adatoms at RT preventing such a charge transfer explains the absence 
of the new band at RT. 

PACS numbers: 73.21.-b,73.22.-f,71.10.Pm 



I. INTRODUCTION 

Recent spur on the study of graphene has been mo- 
tivated mostly by the unusual nature of its effectively 
massless charge carriers so-called Dirac fermions lead- 
ing to some exotic electronic properties, its potential 
appHcation for future nanometer-scale devices, on the 
other hand, has also been under extensive exploration. 
Graphene is a flat atomic single layer of graphite having 
a hexagonal crystalline symmetry with two atoms per 
unit cell. Its low energy band structure shows the almost 
linear tt bands with a vanishing density of states at the 
Dirac point of charge neutrality where the valence band 
(tt band) and conduction band (tt*) meet at the K point 
of the hexagonal Brillouin zone. With increasing reports 
on the electronic properties of graphene layers includ- 
ing the anomalous quantum Hall effects [2- 4], the possi- 
bility of utilizing such features for atomic-scale devices 
seems to increase as seen also in the bilayer graphene 
[14]. Ohta et al. show that the adsorption of Potassium 
on bilayer graphene modifies the potential symmetry in 
the graphene layers resulting in the change of energy gap 
at the K point. This provides a way of controlling the 
conductivity of the charge carriers for the switching capa- 
bility in the atomic scale. Moreover recent calculations 
also provide a way to control the electrical properties 
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of graphene by utilizing the chiral nature of the Dirac 
fermions of graphene. [5] Park et al., show that the Dirac 
fermions of graphene can be quite sensitive to an ex- 
ternal potential resulting, for example, the anisotropic 
renormalization of group velocity at the K point driven 
essentially by the chiral character of charge carriers of 
graphene. This feature becomes more prominent when 
the potential becomes periodic with a periodicity L com- 
parable to that of the graphene lattice constant ao ~1.42 
A. In addition, the life-time of Dirac fermions becomes 
significantly reduced when L<Cao by enhancing back- 
scatterings from the potential between states with mo- 
mentums of -\-k and —k [sl, 0- No experimental study on 
the stability of graphene bands against a foreign pertur- 
bation, however, has been made especially for a potential 
with L~ao despite several earlier angle-resolved photoe- 
mission spectroscopy (ARPES) studies (sl-IToj. 

We report here an example of such a study that the 
presence of an external potential on graphene signifi- 
cantly affects the nature of Dirac fermions not only by 
the periodicity of the potential but also more remarkably 
by the degree of order of the potential especially when 
L~ao. We have investigated the behavior of the tt band 
along the T — K direction as we adsorb Na on graphene 
at two different temperatures, 85 K and room tempera- 
ture (RT). We find remarkable changes in the tt band as a 
function of Na dose at both temperatures. We also notice 
the appearance of a parabolic band centered at k=1.15 
A~^ only at 85K. Furthermore the tt band exhibits quite 
different behavior with Na on a bi-layer graphene from 
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that of a single layer graphene. We mainly discuss such 
changes in the tt band of graphene in terms of the varia- 
tion of the induced potential by Na adsorption. We also 
discuss the origin and characteristics of the new parabolic 
band from our band calculation. 



II. METHODS 

We have obtained our ARPES data by using the 
ARPES chamber at the beamline 3A2 of the Pohang 
Accelerator Laboratory in Korea using synchrotron pho- 
tons of energy 34 eV. The chamber was equipped with a 
hemispherical Scienta R4000 electron analyzer which pro- 
vides an overall energy resolution of 113 meV. All mea- 
surements were made with the chamber under lxlO~^^ 
Torr. An n-type 6H-SiC(0001) substrate was pre-treated 
by annealing at 900 °C under the Si flux generated from 
a resistively heated Si wafer at 1150 °C to obtain the 
>/3 X >/3 phase. We then continued to anneal the sub- 
strate further to obtain a carbon rich surface with the 
6a/3 X 6a/3 phase as reported earlier. [? ] Single layer and 
bilayer graphenes were self-assembled when this precur- 
sor phase was annealed for an extended time period. We 
have used a commercial (SAES) getter source to deposit 
Na on the graphene surface and measured work function 
change to estimate Na coverage. 

In order to interpret our band data, we have also 
performed ab initio pseudopotential total energy calcu- 
lations with a plane- wave basis set. [21] The exchange- 
correlation of electrons was treated within the general- 
ized gradient approximation (GGA) as implemented by 
Perdew-Berke-Enzelhof.[22j The cutoff energy for the ex- 
pansion of wave-functions and potentials in the plane- 
wave basis was chosen to be 400 eV, and the Bril- 
louin zone sampling was done with the Monkhorst-Pack 
special k-point method with a grid of 5x5x1 for the 
a/S X buffer layer phase. We used the projector aug- 
mented wave pseudopotentials as provided by the soft- 
ware package (Vienna Ab-initio Simulation Package). [23] 
The atomic relaxation was carried out until the Helmann- 
Feynman forces were less than 0.02 eV/A. The vacuum 
layer in the supercell used in our calculations is set to be 
8 A, which is enough to minimize the artificial inter layer 
interaction. 



III. RESULTS AND DISCUSSIONS 

The phase of single layer graphene formed on 
SiC(OOOl) surface was first indicated by its characteristic 
LEED spots as seen in Fig. [TJa). While weakening the 
spots from the underlying a/3 x a/3 phase (arrow 1), the 
LEED spots unique to the graphene 1x1 phase (arrow 
2) appear with growing intensity as annealing continued. 
A primary electron energy of 98 eV is used to monitor 
the changes in the LEED pattern during the prepara- 
tion process. The electronic structure of the graphene 




FIG. 1: (Color online) (a) A typical LEED pattern from a 
single layer graphene formed on SiC substrate measured with 
a primary beam energy of 98 eV. The characteristic spots (ar- 
rows 1 and 2) from the single layer graphene coexist with the 
spots (arrow 1) from the VS x y/3 phase, (b) The formation 
of a single layer graphene is finally confirmed by observing the 
unsplit linear tt band of graphene in the ARPES intensity map 
near the K point. The band was measured with synchrotron 
photons of energy 34 eV along the F — K direction. 



thus prepared has been measured using ARPES along 
the T — K direction of graphene up to the second Bril- 
louin zone of the V3 x ^/3 phase [see Fig. Efb)]. Photo- 
electron intensity map near the zone boundary K shows 
the well known graphene tt band [Fig. [Hb)] an indica- 
tive of a single layer graphene. This linear dispersion of 
the TT band suggests that the charge carriers in graphene 
can be considered as massless Dirac fermions. The single 
layer graphene suggested by LEED, however, has been 
finally confirmed by observing the unsplit linear tt band 
along the T — K direction as shown in Fig. [Hb) before 
developing into the split tt band with further annealing 
due to the formation of multiple graphene layeres.[14| 
Here we observe only one band rather than the two form- 
ing a conical shap e tt band due to the chiral property of 
graphene. [1^ Il3| From the slope of the band we esti- 
mate the Fermi velocity Vi7=l.l±0.074x 10^ m/s in good 
agreement with earlier value of 10^ m/s from Shubnikov- 
de Haas oscillations pL| where the effective carrier mass of 
graphene was estimated to be about 0.02 ^ 0.07 mo with 
mo being the free electron mass.[l] 

In Fig. [21 we show progressive changes of the tt band 
with increasing function of Na dose on a single layer 
(a-e) and a bi-layer (f-j) graphene at temperature 85 
K. The amount of Na atoms adsorbed is indicated by 
work function change (A0) induced by Na adsorption. 
The line- width of the tt band in (f) appears to be much 
broader than that in (a) of the single layer graphene. 
Although our resolution does not resolve the broadened 
band into two or more split bands, the measured line- 
width of 0.08=b0.02 A~^ indicates the tt band from a bi- 
layer graphene. [l^ From (a-e), one immediately notices 
that the adsorption of Na on the single layer graphene 
quickly deteriorates the tt band by significantly weaken- 
ing its spectral intensity as well as broadening the line- 
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FIG. 2: (Color online) (a)-(e) The progressive change of the band structure of a single layer graphene along the T — K direction 
as a function of Na dose at 85 K. With increasing Na dose, the linear tt band quickly deteriorates both in intensity and in band 
width while a new parabolic band develops, (f)-(j) Corresponding change for a bi- layer graphene at 85 K. 



width, and concomitantly accompanies a new parabolic 
band near Fermi energy. The deterioration of the tt 
band becomes worse while the parabolic band becomes 
stronger with increasing Na coverage. The correspond- 
ing changes induced by Na adsorption on the bi-layer 
graphene (f-j) are quite different from those on the single 
layer graphene; The tt band remains strong in intensity 
while slightly broadening with increasing Na dose. The 
parabolic band is still there but becomes much weaker 
compared to the one on the single layer graphene at the 
same dose of Na. 

We first investigate the nature of the parabolic band. 
Figure [3] shows this new band at A(/)=0.30 eV. Red cir- 
cles in Fig. [3fa) are the maxima in spectral intensity ob- 
tained by fitting the energy distribution curves (EDCs) 
with a Lorentzian function for 0.98 <k<1.32 
[see Fig. ISl^c)]. The dispersion of the new band esti- 
mated from the momemtum distribution curve at Fermi 
energy in Fig. [3jb) appears to be parabolic with a min- 
imum at k ~1.15 for E=0.19 eV below the Fermi 
level. We obtain charge carrier mass m=0.58±0.07 mo 
and i'i?=3.4±0.41xl0^ m/s, from the parabolic disper- 
sion of the new band. The effective mass of charge carri- 
ers of this band is 8.3^29 times heavier than that of the 
Dirac electrons of the clean graphene without Na. 

Interestingly, the tt band quickly fades away with 
the development of the new band as Na dose increases 
on the single layer graphene, which is more apparent 
in Fig. m^a) where the momentum distribution curves 
(MDCs) change with Na dose. The tt band becomes al- 
most disappeared into the broad background, whereas 



the new band enhances its spectral intensity showing 
the two peaks marked by arrows of the parabolic band 
centered at around 1.15 A . The tt band, however, be- 
haves quite differently with Na for the bi-layer graphene 
mostly maintaining its intensity as seen in Fig. HJ^b). 
The rapid quenching of the tt band of the single layer 
graphene reflects a destructive effect associated with the 
Na adatoms on the top graphene layer. The sustaining 
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FIG. 3: (Color online) (a) Photoemission intensity map of 
a new band from a single layer graphene along the T — K 
direction at Na dose of A(/)=— 0.30. (b) The new band appears 
to be a parabolic shape with a minimum at^l.l5A~^ , which 
is apparent in the MDC at the Fermi level and also in (c) 
where EDCs are shown between ki=0.98 and A:/ = 1.32 
A~^. The red circles denote the maximua in intensity. 
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FIG. 4: (Color online) The spectral changes of the MDCs with 
increasing Na dose measured in terms of A(f) obtained from 
(a) single layer graphene and (b) bi-layer graphene. The down 
arrows indicate the intensity maxima of the new parabolic 
band. The solid gray curves are Lorentzian fit-curves of the 
TT band of graphene. The vertical dotted line in (b) shows 
the shift of the tt band toward higher binding energy with 
increasing Na dose. 



intensity of the tt band from the bi-layer graphene in (f- 
j) may thus be attributed to the contribution from the 
second layer graphene overcoming the loss from the top 
graphene layer. This idea seems to be supported by the 
change of band width as a function of Na dose shown 
in inset of Fig. Hl^b). The band width decreases initially 
down to the value of a single layer graphene and then be- 
gins to increase with further dose of Na. We will discuss 
later the possible causes of such changes of the tt band 
in Fig. El 

We also performed ab initio pseudopotential total 
energy calculations with plane- wave basis set. [21] The 
exchange-correlation of electrons was treated within the 
generalized gradient approximation (GGA) as imple- 
mented by Perdew-Berke-Enzelhof.[22j The cutoff energy 
for the expansion of wave-functions and potentials in the 
plane-wave basis was chosen to be 400 eV, and the Bril- 
louin zone sampling was done with the Monkhorst-Pack 
special k-point method with a grid of 5x5x1 for the 
v^S X buffer layer phase. We used the projector aug- 
mented wave pseudopotentials as provided by the soft- 
ware package (Vienna Ab-initio Simulation Package). [23] 
The atomic relaxation was carried out until the Helmann- 
Feynman forces were less than 0.02 eV/A. The vacuum 
layer in the supercell used in our calculations is set to be 
8 A, which is enough to minimize the artificial interlayer 
interaction. 

Figure [5] shows the plausible binding site of Na atoms 
on the ^/3 x buffer layer in (a) and the Brillouin zones 
of the 1 X 1 of graphite together with those of the buffer 
layer in (b). With the Na coverage of 0.63 ML, the ad- 



FIG. 5: (Color online) (a) Side (left) and top (right) view of 
the atomic arrangement for the >/3 x buffer layer/6H-SiC 
assuming Na coverage of 0.63 single layers, (b) The Brillouin 
zone of 1x1 clean graphene (red hexagon) and that of the 
VS X a/3 buffer layer (blue hexagon). The second Brillouin 
zone of the a/3 x y/3 phase shifted by 1.20 is also shown 
by pale blue hexagon. Calculated band structures (c) for the 
a/3 X a/3 buffer layer/6H-SiC system and (d) for the single 
layer graphene/buffer layer/6H-SiC system. Note that the tt 
bands of the buffer layer near Fermi energy in the vicinity of 
r point appeared to be partially occupied in (c) for the buffer 
layer while they are empty in (d) for the graphene layer. 



sorption energy Ead calculated for the most favored site 
turns out to be 0.86 eV, compared to the next favored 
hollow site of 0.72 eV. All hollow sites of graphene ap- 
pear to have £^^=0.73 eV. We note that the SiC surface 
is not completely covered with either single- or bi-layer 
graphenes.[17j We find the most stable bonding config- 
uration with a maximum adsorption energy £^^=1-03 
eV when Na atoms form a single crystalline layer cor- 
responding to the (110) plane of BCC Na as also seen in 
the scanning tunneling microscopy measurements of Na 
islands on graphite. |15| 

In Figure O we also present our band structures calcu- 
lated with Na-added on the buffer layer in (c) and on the 
graphene layer in (d). One immediately finds the three 
partially filled bands near the F point in (c) originating 
the upper two from the tt electrons in the buffer layer and 
the lower one from the Na adatoms. Since the buffer- layer 
TT bands are completely empty without Na adatoms [2^, 
the TT bands of the buffer layer are rigidly shifted down 
by charge transfer from Na adatoms to the tt bands. The 
nearly empty upper tt bands in (d) then indicate that the 
charge transfer from Na occurs only to the top graphene 
layer with the underlying buffer layer almost intact. As 
depicted in Figure [5] (b) since the separation between the 
first (Fo) and second (Fi) Brillouin zones of the buffer 
layer is Ak=1.20 A~^ much shorter than Ak=2.10 A~^ 
of the Na (110) plane, the new parabolic band observed in 
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Figure [2] must be the tt bands of the buffer layer partiahy 
filled by electrons from Na adatoms. The effective masses 
of the two buffer-layer tt bands are calculated to be 0.63 
mo and 1.35 mo, respectively. We note that the effec- 
tive mass of the upper one is quite close to the observed 
value of 0.58±0.07 mo. The relatively weak parabolic 
bands seen from the bilayer graphene in Figure [5] (f-j) 
may be caused by the reduced area of the exposed buffer 
layer compared to that of the single layer graphene. 

We now consider why the linear tt band of graphene 
deteriorates with Na adsorption. One may think of two 
plausible causes ; 1) due to the random scattering of 
photoelectrons from the disordered Na adatoms or 2) 
due to the reduced mean-free path of charge carriers by 
enhanced back-scatterings of charge carriers from a Na- 
induced periodic potential. The first cause, however, can 
be easily ruled out by our data in Fig. [6] showing the cor- 
responding behavior of the tt band of graphene at RT. 
Here one immediately notices neither a new parabolic 
band nor the deterioration of the tt band with increasing 
Na. Instead, the tt band becomes stronger in intensity 
and shifts toward the higher binding energy by about 0.5 
eV. Interestingly we find that the strong tt band at RT 
begins to deteriorate upon cooHng the sample eventually 
reproducing the low-temperature features in Fig. [2l The 
change in the tt band of graphene with Na appears to 
be reversible in temperature. We thus attribute the sec- 
ond cause to the Na-induced weakening of the tt band 
of graphene assuming that Na adatoms form crystalline 
islands at low temperature giving rise to a certain form 
of a periodic potential. 

If Na crystalizes at 85 K in the form of a layer cor- 
responding to the (110) plane of BCC Na crystal as the 
Na islands formed on graphite, [15] the periodicity L of 
the Na-induced potential may be similar to ao of the 
graphene lattice constant. Such a potential may severely 
modify the band structure in a way to reduce the mean- 
free path of charge carriers within graphene by enhanced 
back-scatterings of charge carriers as noted earlier. [5|- 
0, [13] The formation of such crystalline Na islands at 
85 K is indeed highly likely considering the significantly 
reduced hopping rate between the hollow binding sites 
(z/ = z^oe"^^/^^-^, where is an attempt frequency and 
kB is the Boltzmann constant), which is about 3.65x10^ 
times smaller than that at RT. Here we have used the 
diffusion barrier of 0.13 eV between the most stable hol- 
low sites for Na on graphene. The stronger and shifted 
TT band at RT suggests that the tt band is shifted down 
rigidly towards the higher binding energy side revealing 
the highly occupied band near Fermi level as also ob- 
served when K is absorbed. [25] On the other hand, the 
increased hopping rate of Na at RT drives delocalization 
of charge carriers [23 and reduces the donation of charges 
to the buffer layer leaving the tt band of the buffer layer 
empty as we find Fig. [6l 
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FIG. 6: (Color online) (a)-(d) The corresponding changes as 
in (a)-(e) of Fig. 2 with the sample at RT. The changes are 
quite different from those in Fig. 2 showing the tt band rigidly 
shifted toward the higher binding energy with enhanced spec- 
tral intensity. One also notes the absence of a new parabolic 
band near Fermi energy in the vicinity of F. 



IV. CONCLUSION 

By adsorbing Na on graphene at two different temper- 
atures (RT and 85 K) we observe that the band structure 
of graphene formed on SiC substrate becomes quite sen- 
sitively modified by the Na-induced periodic potential es- 
pecially when the periodicity of the potential is compara- 
ble to the lattice constant of graphene. The frozen tt band 
of graphene quickly weakens in intensity and broadens in 
line- width with Na added at low temperature while it be- 
comes stronger but is solidly shifted towards the higher 
binding energy with Na at RT. A new parabolic band 
observed only at 85 K by adding Na is found to be the tt 
band of the underlying buffer layer partially filled by the 
charges from Na. Such an extremely sensitive nature of 
Dirac fermions in graphene to external impurities as well 
as temperature demands a more careful study of trans- 
port properties of graphene for future practical applica- 
tions in developing graphene based electronic devices. 
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